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Abstract. We construct explicit examples of one-dimensional driven diffusive sys- 
tems for two and three species of interacting particles, defined by asymmetric 
dynamical rules which do not obey detailed balance, but whose nonequilibrium 
' O ' stationary-state measure coincides with a prescribed equilibrium Gibbs measure. For 

simplicity, the measures considered in this construction only involve nearest-neighbor 
^ ' interactions. For two species, the dynamics thus obtained generically has five free 

parameters, and does not obey pairwise balance in general. The latter property 
pv^ I is satisfied only by the totally asymmetric dynamics and the partially asymmetric 

K^ ' dynamics with uniform bias, i.e., the cases originally considered by Katz, Lebowitz, and 

2]H ■ Spohn. For three species of interacting particles, with nearest-neighbor interactions 

*^ ' between particles of the same species, the totally asymmetric dynamics thus obtained 

^^ ■ has two free parameters, and obeys pairwise balance. These models are put in 

^^ , perspective with other examples of driven diffusive systems. The emerging picture is 

^— s ' that asymmetric (nonequilibrium) stochastic dynamics leading to a given stationary- 

state measure are far more constrained (in terms of numbers of free parameters) than 
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1. Introduction 

Driven diffusive systems [1, 2, 3, 4] are defined by stochastic dynamical rules that 
incorporate the effect of an external drive, and therefore do not obey detailed balance, 
which makes their nonequilibrium stationary states difficult to study in general. The 
first step of evaluating the stationary-state measure, i.e., the probability Pst(C) of any 
given configuration C, is already a non-trivial task. In view of the lack of a general 
theory of nonequilibrium stationary states, one has to rely on the investigation of specific 
models, for which the stationary-state measure is analytically tractable. Examples of 
such models are the simple exclusion processes, the zero-range process, and the one- 
dimensional Katz-Lebowitz-Spohn (KLS) model. 

For the simple exclusion process, particles obeying an exclusion constraint perform 
random walks on a lattice, with either symmetric or biased moves [3, 4, 5]. On a 
ring, i.e., with periodic boundary conditions, the stationary-state measure is uniform, 
irrespectively of the bias: all the configurations are equally probable. The same model 
with open boundaries has a stationary-state measure which can be described in terms 
of a matrix product Ansatz. The zero-range process (ZRP) [6], in the presence of 
a bias, belongs to the class of driven stochastic processes with multiple occupancies. 
Its stationary state has a product measure, again irrespectively of the bias. The 
occupation numbers of the sites are independent random quantities with a common 
distribution, up to the conservation of the total number of particles. The existence 
of unbounded occupancies however opens up the possibility of having a condensation 
transition, irrespective of the geometry of the underlying lattice, and therefore also in 
one dimension [7, 8]. 

The KLS model is a lattice gas model of interacting charged particles subjected to an 
external electric field [9]. It is representative of the class of models with non-equilibrium 
stationary state measures incorporating physical interactions between particles. The 
stationary-state measure of the KLS model is non trivial in two dimensions and above, 
where e.g. the critical temperature depends continuously on the applied field [2]. 
The situation however simplifies in the one-dimensional geometry, where the model 
is equivalent to a chain of classical Ising spins s„ = ±L In this case, there exists a class 
of biased stochastic dynamics, for which the nonequilibrium stationary-state measure 
is the canonical finite-temperature Gibbs measure, where the probability Pst(C) of the 
configuration C = {si, . . . , Sn} is given by the Boltzmann formula (with ksT = 1) 

Pst(C) = ^ exp{-n{C)) (LI) 

associated with the Ising Hamiltonian with nearest-neighbor interactions 

7i = -J^s„s„+i. (1.2) 

n 

This very stochastic model with antiferromagnetic interactions (J < 0) was subsequently 
rediscovered in the context of polymer crystallization [10]. 

At this point it is natural to question the generality of the result found in [9] for the 
case of the Ising chain. In the present work we show that for systems with three species 
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of particles there also exist asymmetric stochastic dynamics which do not obey detailed 
balance, but whose stationary-state measure is the Gibbs measure corresponding to a 
simple Hamiltonian. In particular this stationary-state measure is independent of the 
bias. In our construction we restrict the choice of measures to those involving only 
nearest-neighbor interactions. We first revisit in Section 2 the case of two species of 
interacting particles, thus generalizing the study of [9] to a wider class of dynamics. 
The emphasis is put on the role of various symmetries, and especially on the number 
of free parameters left by imposing them. We then consider, in Section 3, the entirely 
novel situation of three species of interacting particles. Section 4 contains a discussion, 
where our results are put in perspective with yet other examples. 

Let us finally give a brief reminder of the concepts of detailed balance [11] and 
pairwise balance [12], which will be used throughout this work. Consider a finite set 
of configurations C, and a Markovian dynamics in continuous time, defined by the 
transition rates W{C -^ C). The master equation for the time-dependent probability 
P(C,t) reads 

^^^^ = E(W^(C' - C)P{C',t) - W{C ^ C')P{C,t)). (1.3) 

at ^1 

The stationary probability Pst(C) therefore obeys the equation 

Y.{W{C' ^ C)Pst(C') - W{C ^ C')Pst(C)) = 0. (1.4) 

c 

In the particular case where the Markov process is reversible, the dynamics 

brings the system to an equilibrium state. Reversibility requires the detailed balance 

property [11], that is the absence of probability flux between any pair of configurations C 

and C at stationarity: 

W{C ^ C')Pst(C) = W{C' ^ C)Pst(C'). (1.5) 

Equation (1.5) clearly implies (1.4). 

A weaker property, referred to as pairwise balance [12], is adapted to the situation 
of driven diffusive systems, where the presence of a preferred direction of motion, i.e., a 
bias, precludes the property of detailed balance. Pairwise balance is defined as follows: 
for every pair of configurations C and C such that W{C ^' C) j^ 0, there exists a third 
configuration C" such that 

WiC ^ C')Pst(C) = W{C" ^ C)Pst(C"). (1.6) 

The moves C ^ C and C" -^ C are said to be conjugate to each other. Equation (1.6) 
also implies (1.4), since terms corresponding to pairs of conjugate moves cancel each 
other. Several of the examples of dynamics constructed in this paper obey P-related 
pairwise balance, where pairs of conjugate moves are related to each other by parity. 
The first explicit example is given below (2.20). An important general consequence of 
P-related pairwise balance is shown at the end of Section 3. If the totally asymmetric 
dynamics obeys P-related pairwise balance, the partially asymmetric dynamics with 
uniform bias p also obeys pairwise balance, and it has the same stationary-state measure 
as the totally asymmetric one, irrespective of the value of the bias. 
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2. Two species 

Consider a ring made of A^ sites. Each site is occupied by a particle, which can be either 
of type A (positively charged) or of type B (negatively charged). We represent the 
species of particle at site n by an Ising spin s„ = ±1 equal to the charge of the particle. 
Table 1 also gives the corresponding indicator variables. For instance, (l + s„)/2 is equal 
to 1 if the particle at site n is of type A, and to else. 



Table 1. Spin (charge) and indicator variables associated with each particle species 
in the case of two species. 



Species at site n 


Spin (charge) 


Indicator variable 


A 
B 


Sn = -1 


(l + Sn)/2 

(l-s„)/2 



Our goal is to construct nonequilibrium dynamics such that the stationary-state 
measure is the Gibbs measure given by (1.1) associated with the nearest-neighbor 
Hamiltonian (1.2). We consider the asymmetric exchange (Kawasaki) dynamics. 
Consistently with the form of the Hamiltonian (1.2), the rates depend on the two 
neighbors of the pair to be exchanged, according to Table 2. The dynamics so defined 
conserves the numbers Na and Nb of particles of each species, with Na + Nb = N. We 
first determine the most general dynamics of this form in Section 2.1, and then discuss 
the interplay between various possible symmetries in Section 2.2. 



Table 2. List of moves in the general exchange dynamics for two species of interacting 
particles, notation for the corresponding exchange rates, and energy difference ATi. 
involved in the moves, where the Hamiltonian Ti. is defined in (1.2). 



Move 


Rate 


AH 


Move 


Rate 


AH 


AABA- 


^ ABAA 


waa 





ABAA- 


^ AABA 


xaa 





AABB - 


^ ABAB 


wab 


4J 


ABAB - 


^ AABB 


xab 


-4J 


BABA- 


-^BBAA 


wba 


-4J 


BBAA- 


-^ BABA 


xba 


AJ 


BABB - 


^ BBAB 


wbb 





BBAB- 


^ BABB 


xbb 






2.1. The general case 

Consider the numbers Naa, ■ ■ ■ , Nbb of oriented pairs of neighbors of each species. These 
numbers obey the sum rules 

Na = Naa + Nab = Naa + Nba, Nb = Nba + Nbb = Nab + Nbb ■ (2.1) 

The sum of the two equations gives twice the same equation, so that the four pair 
numbers obey three independent equations, leaving one single free quantity. It is 
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convenient to take the latter as being the Hamiltonian 7i of (1.2). The pair numbers 
can indeed be expressed as hnear combinations of 7i and of the particle numbers N^ 
and Nb'- 

Naa = \ (3iVA - iVs - HI J) , Nbb = \ (SiVs -Na- Uj J) , 

4 4 (2.2) 

NAB = NBA = \{N + n/J). 

Incidentally, this proves that the Hamiltonian (1.2) is the most general form of a pair 

Hamiltonian for two species of particles with nearest-neighbor interactions. 

Throughout this paper, we make use of an alternative and convenient way of 

automatically taking into account sum rules such as (2.1). This consists in expressing 

all the quantities in terms of the spin variables s„ introduced in Table 1. For instance 

1 A^ 

NA = -T.i^ + Sn) = -{l + {s^)). (2.3) 



n 



Here and in the following, the brackets (...) denote a uniform spatial average for a fixed 
generic configuration C. Recall that all the sites are equivalent, because of translational 
invariance. The pair numbers and the Hamiltonian read 

A^ N 

Naa = j{1 + '2 (si) + (S1S2)), iVfiB = J (1 - 2 (si) + (S1S2)), 

(2.4) 
Nab = iV^A = j (1 - (S1S2)), n = -NJ (siSa) . 

Equations such as (2.1) and (2.2) are then automatically satisfied. 

Consider now the fate of a generic configuration C. The total exit rate Wont{C) 
from C to any other configuration C per unit time can be read off from Table 2: 

Wont{C) = waaNaaba + wabNaabb + wbaNbaba + wbbNbabb 

+ xaaNabaa + xabNabab + xbaNbbaa + xbbNbbab- (2.5) 

An analogous expression can be derived for the total entrance rate Win{C) from any other 
configuration C to C. Using again Table 2, as well as (1.1) to express the stationary-state 
weight Pst(C) as 

Pst(C') = Pst(C)exp(A7^), (2.6) 

in terms of Pst(C) and of the energy difference 

AH = n{C) - n{C') (2.7) 

involved in the move, we obtain 

Wi^{C) = waaNabaa + ^^'^wabNabab + ^'^^'^wba^bbaa + wbbNbbab 

+ xaaNaaba + g'^'^xabNaabb + g'^'^xbaNbaba + xbbNbabb- (2.8) 

In the stationary state we have 

Wont{C)-WUC) = (2.9) 
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for every configuration C. In order to determine the number of independent conditions 
on the rates imposed by this equation, it is convenient to rewrite (2.5) and (2.8) in terms 
of spin correlations, i.e., spatial averages of products of spin variables, denoted as (...), 
along the lines of (2.3). With these notations, we obtain 

Wout(C) - W-UC) = ^ {(e-^-^ - 1)((S1S2S3S4) + l)Rl+ (S1(S3 - S2)S4) i?2 

(si(3s2 - 2s3 + s,)) + e-^' (si(s2 - 2s3 - S4)) 1 i?i j, (2.10) 



+ 

where i?i and R2 stand for the following linear combinations of the rates: 

Ri = e'^'^{wAB + xba) - wba - xab, 

R2 = (e^-^ + l)iwAB - xba) + (e"^-^ + 1)(wba - xab) 

-2{waa + wbb-xaa-xbb)- (2.11) 

The condition (2.9) therefore gives two linear relations, 

R, = R2 = 0, (2.12) 

between the eight exchange rates defining the general asymmetric dynamics. Let us 
choose the time unit by setting 

waa + wbb + xaa + xbb = '^- (2.13) 

The most general asymmetric dynamics for two species of interacting particles 
such that the stationary-state measure is given by (1.1), (1-2) therefore has five free 
parameters. An explicit parametrization of the rates is given in Table 3. The dynamics 
thus obtained does not obey pairwise balance in general. 

The parametrization of the solutions to (2.12) and (2.13) given in Table 3 has 
been carefully chosen in such a way that the various symmetries to be described 
below correspond to the simple constraints (2.17), (2.21), (2.25), (2.27) in terms of 
the parameters 6 and Ei, . . . ,£4. The parameters S and ei, . . . ,£4 all lie in the range 
[—1, +1], and are such that the combination 

A = (l + ^)^i + (l-^)^^ (2.14) 

£3 + ^4 
is positive. 

2.2. The interplay between various symmetries 

The number of free parameters of the dynamics thus obtained is decreased if various 
kinds of symmetries are imposed onto the dynamics. 

• Symmetric (P -invariant) dynamics. Consider a symmetric dynamics, invariant under 
the spatial parity P which reverses the orientation of the ring (i.e., interchanges left and 
right). This symmetry property reads 

xij = wji (2.15) 
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Table 3. Explicit parametrization of the rates of the most general asymnietric 
dynamics for two species with Gibbs stationary-state measure (1-1), (1-2). The 
notation A is defined in (2.14). 



Rate 


expression 


Rate 


expression 


waa 
wab 
wba 
wbb 


4 

(1+^3)A 

2{e^^ + 1) 

(l+£4)Ae^-^ 
2{e^' + 1) 

4 


xaa 
xab 
xba 
xbb 


4 
(l-e4)Ae^-^ 


2{e^' + 1) 

(1-£3)A 

2(e^-^ + 1) 

{l-S2){l-6) 
4 



for all values of the indices I, J = A, B. The stationary state thus obtained is an 
equilibrium state. The first equation of (2.12), 

Wba = e^'^WAB, (2.16) 

expresses detailed balance. Equation (2.15) amounts to setting 

£, = (2 = 1,..., 4) (2.17) 

in Table 3. The symmetric (equilibrium) dynamics therefore has two free parameters: S 
and A, in the ranges — 1 < 5 < 1 and A > 0. The expression (2.14) for A indeed becomes 
indeterminate in the limit where all the Ei go simultaneously to zero. The rates read 



Waa = xaa 



Wba = Xab 



1 + 6 



A 



Ae 



4J 



2(e^^ + 1] 



Wab = xba 



Wbb = Xbb 



2(e'' + 11 



1-S 



(2.18) 



• Totally asymmetric dynamics. Consider a dynamics driven by an infinitely strong 
electric field, so that the positively (resp. negatively) charged A particles (resp. B 
particles) hop exclusively to the right (resp. to the left). Therefore 

xij = (2.19) 

for all values of the indices I, J = A, B. Equation (2.12) becomes 

Wba = e^'^WAB- (2.20) 

This equation coincides with (2.16). It expresses P -related pairwise balance: conjugate 
moves are related to each other by parity P, i.e., the first and the fourth move of the left 
column of Table 2 are their own conjugates, whereas the second and the third moves 
are conjugate to each other. Equation (2.19) amounts to setting 

5, = 1 (z=l,...,4) (2.21) 
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in Table 3, so that A = 1. The totally asymmetric dynamics with stationary-state 
measure (1.1), (1-2) therefore has one free parameter: 6, in the range —1 < 6 < 1. The 
rates read 

1 + S 1 

WAA = ^— , WAB = ^iJ— T' 

^ ^ +1 (2.22) 



^BA = AT . ^ , WbB 



1-5 



e^-^ + r "^^ 2 

• Partially asymmetric dynamics with a uniform bias. This is the most general case 
originally considered by KLS [9]. Consider a dynamics driven by a finite electric field, 
so that the positively (resp. negatively) charged A particles (resp. B particles) hop 
preferentially to the right (resp. to the left). Let 

P = -^, q = -Y- (2-23) 

be the a priori probabilities of respectively hopping to the right and to the left, where 
< e < 1 provides a measure of the applied electric field. This translates into the 
following uniform bias condition: 

^ = 1^ (2.24) 

wji l + e 

for all values of /, J = A, B. This situation interpolates between the symmetric case 

{p = 1/2, e = 0) and the totally asymmetric one {p = 1, e = 1). Equation (2.24) 

amounts to setting 

ei = e (2 = 1,..., 4) (2.25) 

in Table 3, so that again A = 1. As a consequence, there is a two-parameter family of 
dynamics with uniform bias and stationary-state measure (1.1), (1.2), parametrized by e 
and S. We thus recover the original KLS model [9]. The fact that the stationary-state 
weights are independent of the bias is actually a general property of dynamics obeying 
P-related pairwise balance (see Section 3.2). 

• CP -invariance. The CP operation is the product of C and P, where the charge 
conjugation C changes the charge of the particles to its opposite (i.e., interchanges A 
and B particles), whereas the spatial parity P changes the orientation of the ring (i.e., 
interchanges left and right). In physical terms, in the stationary state of a CP-invariant 
dynamics, the current due to a positively charged particle and to a negatively charged 
particle are equal. Requiring CP-invariance yields the two conditions 

waa = wbb, xaa = xbb, (2.26) 

which amount to setting 

5 = 0, £1 = 52 = e (2.27) 

in Table 3. The most general CP-invariant dynamics has therefore three free parame- 
ters: e, £3, and 64. It does not obey pairwise balance in general. 

CP-invariance can be combined with any of the above symmetries: 
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•k The CP-invariant symmetric dynamics corresponds to 6 = Ei = 0. It has a single free 
parameter: A. The rates read 



waa = xaa = T, 



A 



wab = xba 



2(e^^ + 11 



Ae' 



4J 



(2.28) 



wba = xab 



2(e*' + 1] 



WbB = XbB = -7- 



-k The partially asymmetric CP-invariant dynamics with uniform bias has one single free 
parameter: e. 

-k The totally asymmetric CP-invariant dynamics is the most constrained of all the 
dynamics: it has no free parameter at all. The rates 

1 1 e^-^ 

waa = wbb = ^, ^^^ = ^iJ^' ^^^ = ^17^ (2.29) 

only depend on the energy difference A7i involved in the exchange moves. They coincide 
with those of the heat-bath rule [13, 14]: 

1 



w(An) = —— = -(l- tanh ^) . 

^ ' exp{An) + 1 2 V 2 / 



(2.30) 



The above discussion is summarized in Table 4, giving the number of free parameters 
for every symmetry class of dynamics, both without and with imposing CP-invariance. 



Table 4. List of the symmetry classes of dynamics for two species of particles with 
Gibbs stationary-state measure (1.1), (1-2), with balance property: detailed balance 
(DB) or pairwise balance (PB), and number of free parameters, both without and with 
CP-invariance. 



Class of dynamics 


Balance property 


Without CP 


With CP 


General 

Symmetric (equilibrium) 

Totally asymmetric 

Partially asymmetric (uniform bias) 


none 
DB 
PB 
PB 


5 
2 
1 
2 


3 

1 

1 



3. Three species 

Consider again a finite ring of A^ sites. Each site is now occupied by a particle which can 
be either of type A (positively charged), of type B (negatively charged), or of type C 
(neutral, i.e., with no charge). We again represent the species of particle at site n by a 
spin Sn = 0, ±1 equal to the charge of the particle, as shown in Table 5. 

We consider Gibbs measures corresponding to the most general (ferromagnetic or 
antiferromagnetic) Hamiltonian involving pairs of identical nearest neighbors: 

n = -2{JaNaa + JbNbb + JcNcc), (3.1) 



Nonequilibrium stationary states with Gibbs measure for two or three species 



10 



Table 5. Spin (charge) and indicator variables associated with each particle species 
in the case of three species. 



Species at site n 


Spin (charge) 


Indicator variable 


A 
B 
C 




Sn{Sn+ l)/2 
SniSn - l)/2 



where the couphng constants Ja, Jb, and Jc can take both signs. The factor 2 
is introduced for consistency with the case of two species. Using the spin variables 
Sn = 0, ±1 defined in Table 5, the Hamiltonian (3.1) can be rewritten as 



H = Eq — ^22 [("^^ + Jb + 4:Jc)SnSn+i + (Ja — Jb) (Sn + Sn+i) + Ja + Jb 



Jfi^n+li (<J.2j 



where Eq = 2{N — 2Nc)Jc is a constant. This is a generalized Blume-Emery-Grifiiths 
spin-1 Hamiltonian [15]. 

We again address the question of the existence of nonequilibrium stochastic 
dynamics whose stationary-state measure is the measure (1.1) associated with the 
Hamiltonian (3.1). The results obtained in Section 2 for two species of particles suggest 
that the case of totally asymmetric exchange dynamics is already of interest. We 
therefore restrict our investigation to this limiting situation for the time being. The 
positively (resp. negatively) charged A particles (resp. B particles) only hop to the right 
(resp. to the left), whereas the neutral C particles can hop in both directions. The 
exchange rates depend on the two neighbors of the pair to be exchanged, according to 
Table 6. The dynamics so defined conserves the numbers Na, Nb, and Nc of particles 
of each species, with Na + Nb + Nc = N. 



3.1. The CP-invariant case 

Motivated by the form of the results of Section 2 on two species, we first consider the 
CP-invariant case, which can be anticipated to be simpler than the generic one. 
As far as statics is concerned, C-invariance implies 

Ja = Jb = J, Jc = Jo- (3.3) 

The Hamiltonian (3.2) becomes the usual Blume-Emery-Grifiiths Hamiltonian [15] 



H = Eq — 22[i'^ ~^ 2 Jo)S'„S'„+i + J 



b„0; 



n^n+1- 



(3.4) 



As far as dynamics is concerned, CP-invariance yields 12 equalities among the 27 
exchange rates: 

waa = wbb, wac = WCB, wbc = WcA, xaa = Vbb, 

xab = Vab, xac = ycB, xba = Vba, xbb = Vaa, (3.5) 

xbc = ycA, xcA = Vbc, Xcb = VAC, Xcc = Vcc- 
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Table 6. List of moves in the totally asymmetric dynamics for three species 
of interacting particles, notation for the corresponding exchange rates, and energy 
difference AH involved in the moves, where the Hamiltonian H is defined in (3.1). 



Move 


Rate 


AH 


Move 


Rate 


AH 


AABA- 


^ ABAA 


waa 





BCBC - 


-^ BBCC 


xbc 


— 2Jb — 2J(j 


AABB- 


-^ ABAB 


wab 


2Ja + 2Jb 


CCBA- 


-^CBCA 


XCA 


2Jc 


AABC - 


^ABAC 


wac 


2Ja 


CCBB- 


^CBCB 


XCB 


2Jb + 2Jc 


BABA- 


^BBAA 


wba 


-2Ja - 2Jb 


CCBC - 


^CBCC 


xcc 





BABB - 


^ BBAB 


Web 





AACA- 


^ ACAA 


Vaa 





BABC - 


^ BBAC 


wbc 


-2Jb 


AACB - 


^ ACAB 


Vab 


2Ja 


CABA- 


^CBAA 


wca 


-2Ja 


AACC - 


^ ACAC 


Vac 


2Ja + 2Jc 


CABB- 


^CBAB 


WCB 


2Jb 


BACA- 


^ BCAA 


Vba 


-2Ja 


CABC - 


^CBAC 


wcc 





BACB - 


^ BCAB 


Vbb 





ACBA- 


^ ABCA 


xaa 





BACC - 


^ BCAC 


Vbc 


2Jc 


ACBB- 


-^ ABCB 


xab 


2Jb 


CACA- 


^CCAA 


ycA 


-2Ja - 2Jc 


ACBC - 


^ ABCC 


xac 


-'iJc 


CACB- 


^CCAB 


ycB 


-2Jc 


BCBA- 


^BBCA 


xba 


-2Jb 


CACC - 


^CCAC 


ycc 





BCBB - 


^ BBCB 


xbb 














The analysis follows the lines of Section 2. The algebra is however far more cumber- 
some, so that intermediate expressions are too lengthy to be reported here. Calculations 
have been worked out with the help of the software MACSYMA. We start from the 
expressions for the total rates WoutiC) and VFin(C) for a generic configuration C, similar 
to (2.5) and (2.8), which can be read off from Table 6. The difference VFout(C) — VFin(C) is 
then recast in terms of products of the spin variables S'„. We thus obtain an expression 
similar to (2.10), involving 42 different correlations of two to eight spin variables. One 
example of a correlation of two variables is {S1S2), whereas there is a unique correlation 
of eight variables: {SfS2S^Sl). Requiring that the coefficients of all these correlations 
vanish, we thus obtain 42 linear (but not independent) relations of the form 

i?i = . . . = i?42 = 0, (3.6) 

where the Ri are linear combinations of the exchange rates, similar to Ri and R2 given 
in (2.11). In the CP-invariant situation under study, equations (3.5) and (3.6) together 
yield 26 independent linear relations among the 27 exchange rates, and therefore leave 
a single free parameter, which can be fixed by choosing the time unit. For consistency 
with the case of two species, we set 

Waa + wbb = 1- (3.7) 

This normalization condition uniquely determines all the exchange rates. 

We have therefore shown that there is a single totally asymmetric CP-inva- 
riant dynamics for three species of interacting particles with stationary-state meas- 



Nonequilibrium stationary states with Gibbs measure for two or three species 



12 



ure (1.1), (3.1). This uniquely determined dynamics can be viewed as a non-trivial 
extension to three species of the result (2.29). The explicit expressions of the exchange 
rates are given in Table 7. This dynamics obeys pairwise balance (see equation (3.11) 
below for the general case). However, at variance with the case of two species, the rates 
are not of the heat-bath form (2.30). 



Table 7. Expressions of the rates of the totaUy asymmetric CP-invariant dynamics 
for three species. The label / stands for any particle species (/ = A, B, C). 



Rate 


expression 


Rate 


expression 


Rate 


expression 


waa = wbb 
wab 

WAC = WcB 


1 
2 

1 


wba 

wbc = WCA 

wcc 


e^-' 


XIA = VBI 
XlB = VAI 

xic = yci 


e'' 


e^-^ + 1 

e''{e^' + 1) 

2{e^' + 1) 


2(e^^ + 1) 
1 


e^^ + 1 

e^^ + l 

2(e^-^ + 1) 


2(e^^ + 1) 

g2J+2Jo 


e^-' + l 


2(e4^ + 1) 



3.2. The general case 

We now turn to the general totally asymmetric dynamics. We view J a, Jb, and Jc 
as three independent coupling constants, and consider the 27 exchange rates entering 
Table 6 as being a priori all different from each other. 

Following the above procedure, and choosing time units according to (3.7), we are 
left after some lengthy algebra with a two-parameter family of dynamics with stationary- 
state measure (1.1), (3.1). An explicit parametrization of the rates is given in Table 8, 
where the parameters a and (3 enter linearly, and with the notation 

f = 2(e2JA+2Js + i)- (2-^) 

The form of the CP-invariant rates given in Table 7 has been helpful in working out 
this parametrization of the general case. 

The two parameters a and (3 run over some domain D, such that all the rates of 
Table 8 are positive. It can be checked that D is an asymmetric quadrilateral, shown 
schematically in Figure 1. The co-ordinates of its vertices read 



"K 



/5k = 0, 



ttL 



Pi 



La 


+ Sb 


La 


-Sb 


2 





ttM 



ftN 



(3u — 0, /5^ 



Lb + Sa 
Lb — Sa 

2 



(3.9) 



L 



B 



Sa 



La — Sb 

where La (resp. Sa) is the largest (resp. the smallest) of the three quantities exp(2J^-f- 
2J(;), exp(— 2J^ + 2J(7), and exp(— 2J4 — 2Jc), and Lb (resp. Sb) is the largest (resp. the 
smallest) of the three quantities exp(2JB+2Jc), exp(— 2Js+2Jc7), and exp(— 27^— 2Jc). 
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Table 8. Explicit paramctrization of rates of the most general totally asymmetric 
dynamics for three species. The notation / is defined in (3.8). 



Rate 


expression 


Rate 


expression 


waa 


(l-a + (l + a)e2'^^+2'^^)/ 


xbc 


(1 + a + /5)e2^«+2'^c/ 


wab 


2/ 


XCA 


iil + a)e'''-+(3)f 


wac 


(l-a + (l + a)e2^s)/ 


XCB 


{l + a + f3)f 


wba 


2g2 Ja +2 Js J 


xcc 


((l + a)e2^s+2Jc+^)/ 


wbb 


(l + a + (l-a)e2-^^+2Js)j 


Vaa 


(1 - a - f3e^J^+^-^c)f 


wbc 


(l + a + (l-a)e2-^^)e2-^sj 


Vab 


{l-a-f3e^^^')f 


WCA 


(l-a + (l + a)e2-^s)e2^-^/ 


Vac 


il-a-(3)f 


WCB 


{l + a + {l-a)e^''^)f 


Vba 


(1 - a - pe^-^^y-^^f 


wcc 


{{l-a)e^^^ + {l + a)e^-^^)f 


Vbb 


{{l-a)e^-^^-f3e^^^)f 


xaa 


{{l + a)e^J^+f3e^'o)f 


Vbc 


iil-a)e'^--P)f 


xab 


il + a + pe^-^^)f 


ycA 


(1 - a - py^^^+^-^cf 


xac 


((l + a)e2-^s+/5)e2^c/ 


ycB 


(^{l-oi)e''JA_i3yjcf 


xba 


(l + a + /5e2^c)e2^^/ 


ycc 


((l_«)e2^A+2Je_^)/ 


xbb 


(l + a + /5e2^s+2-^^)/ 








Figure 1. Typical shape of the quadrilateral domain D of the a-/3 plane such 
that all the rates of Table 8 are positive (see text). 



The general dynamics of Table 8 contains as special cases several of the situations 
considered so far. The CP-invariant situation corresponds to a = /5 = 0. The rates 
of Table 7 are thus recovered, with the notations (3.3). The case of two species is also 
recovered. The rates (2.22) are reproduced, again with the notations (3.3), and with 
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the identification 

e^^-1 
S = ^-j a. (3.10) 

The family of dynamics described in Table 8 obeys P-related pairwise balance. This 
property is expressed by the nine relations 

wba = e^'^-^+^'^swAB, wcA = e'^'^^WAc, wbc = e'^-^'^wcB, 

XBA = e^^^XAB, XAc = e^'^''xcA, XBc:=e2-^^+2-^^a;cB, (3.11) 

Vba = e^'^^VAB, ycA = e^-'-^+^-^^i/Ac, Vcb = e'^-^'^ysc, 

which are identically fulfilled by the rates of Table 8, for any values of the parameters a 
and f3. In other words, the relations (3.11) are built in as a subset of (3.6). 

Pairwise balance has the following consequence. Consider the partially asymmetric 
dynamics with uniform bias p, where all the 'right' moves, i.e, those of Table 6, take 
place with rates equal to p times those given in Table 6, whereas the P-related 'left' 
moves take place with rates equal to q times those given in Table 6, with q = 1 —p. For 
example: 

AABC -^ ABAC with rate pwAc, 

CBAA -^ CAB A with rate qwAc- (3.12) 

For this uniformly biased dynamics, the total entrance and exit rates for a given 
configuration read 

WUC) = PW:^^\C) + qW^iiC), (3.13) 

with self-explanatory notations. The stationarity condition for the totally asymmetric 
dynamics (p = 1) reads H^;"^ (C) = W^u^ (C). On the other hand, P-related pairwise 
balance implies W[^^^'\C) = W^^^^iC) and W^^^\C) = W^f{C). We have therefore 

WUC) = Wi:^''\C) = Wt'\C) = W^out(C) = W:t\C) = W^r (C). (3.14) 

These equations show that the partially asymmetric dynamics with uniform bias p 
has the same stationary-state measure as the totally asymmetric one. This dynamics 
interpolates between the symmetric (equilibrium) case {p = 1/2) and the totally 
asymmetric one {p = 1). The fact that the stationary-state measure is independent 
of the bias p thus appears as a general consequence of P-related pairwise balance. 

4. Discussion 

In this paper we explicitly constructed classes of nonequilibrium dynamics for two 
and three species of interacting particles, i.e., asymmetric stochastic dynamics which 
do not obey detailed balance, but whose nonequilibrium stationary-state measure is a 
prescribed measure. We have chosen to work with finite-temperature canonical Gibbs 
measures associated with spin Hamiltonians with nearest-neighbor interactions. The 
stationary current, as well as many other observables in the stationary state, can 
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therefore be evaluated, at least in principle, by means of the transfer-matrix formalism. 
We have emphasized the role of the various symmetries which can be imposed onto the 
dynamics. 

For two species of interacting particles, a situation first considered by KLS [9], 
stationary-state measures are associated with the usual (ant i) ferromagnetic Hamiltonian 
on the spin-1/2 Ising chain. Our result for the most general dynamics is given by Table 3. 
This dynamics has five free parameters, and does not obey pairwise balance in general. 
Only the cases considered by KLS, namely the totally asymmetric dynamics and the 
partially asymmetric dynamics with a uniform bias, obey P-related pairwise balance, 
where pairs of conjugate moves are related by parity. 

We then turned to the novel situation of three species of interacting particles. 
Stationary-state measures are given by the most general Hamiltonian involving pairs of 
neighboring particles of the same species. This translates into a Blume-Emery-Grifiiths 
spin-1 Hamiltonian. We first restricted the search of dynamics to the totally asymmetric 
case. The most general situation is described by Table 8. This dynamics obeys pairwise 
balance. It can therefore be extended to a partially asymmetric dynamics with uniform 
bias p. The three-parameter family of dynamics thus obtained interpolates between the 
symmetric (equilibrium) case {p = 1/2) and the totally asymmetric one {p = 1). 

The most constrained class of stochastic dynamics we have investigated is the CP- 
invariant totally asymmetric one. For a prescribed stationary-state measure, there is 
indeed a uniquely defined such dynamics, with no free parameter, both for two species 
(see equation (2.29)) and for three species (see Table 7) of interacting particles. 

Throughout this work we have put a strong emphasis on the numbers of free 
parameters in symmetric and asymmetric dynamics leading to a given stationary- state 
measure. Our results suggest the following rule: asymmetric stochastic dynamics 
leading to a given nonequilibrium stationary-state measure are far more constrained 
than symmetric dynamics leading to the same measure as an equilibrium measure. 

To close up, let us demonstrate that the above empirical rule holds in a much 
broader class of stochastic models. To do so, we have chosen to put the results of this 
work in perspective with the following two characteristic examples, which also belong 
to the realm of driven diffusive systems. 

Example 1. 

Our first example is much in the spirit of the present paper. Consider a driven 
diffusive system consisting of K species of non-interacting particles on a ring, denoted 
by / = A, 5, . . ., where K >2 is arbitrary. 

The most general exchange dynamics is defined by the K{K — 1) rates ujj 
corresponding to the moves IJ -^ JI for I ^ J . We look for dynamics such that 
the stationary-state measure is uniform, i.e., all the configurations with given particle 
numbers Nj of each species are equally probable. This is indeed the right concept for a 
Gibbs measure in the absence of interactions, or, equivalently, in the limit of an infinite 
temperature. 
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The condition for having a uniform stationary-state measure reads 

W^out(C) - WUC) =Y.^ij{Nij - Njj) = (4.1) 

ij 

for every configuration C. The number of independent conditions on the rates imposed 

by this equation can be evaluated as follows. There are K^ numbers of oriented 

pairs Njj, which obey the sum rules 

Y.Nij = T. Nji = Ni^ E ^/ = ^- (4-2) 

J J I 

Only {K — 1)^ pair numbers are therefore linearly independent. Equation (4.1) shown 
that each of the {K — 1){K — 2)/2 antisymmetric combinations of these independent pair 
numbers yields one condition. The K{K — 1) rates therefore obey {K — ^){K — 2)/2 
conditions, so that the general asymmetric exchange dynamics for K species of non- 
interacting particles depends on 

AK = liK-l){K + 2) (4.3) 

dimensionful parameters. On the other hand, for the symmetric exchange dynamics 
obeying the detailed balance property uu = uji, the K{K — l)/2 rates are not 
constrained at all. Indeed (4.1) vanishes identically. The general symmetric exchange 
dynamics therefore depends on 

Sk = \K{K - 1) (4.4) 

parameters. One has 

Ak = Sk + K-1. (4.5) 

For two species {K = 2), we have A2 = 2 and S2 = 1- There is no condition on 
the exchange rates, because there exists no antisymmetric combination of pair numbers. 
Equation (2.2) indeed implies Nab = Nba- As a consequence, the stationary-state 
measure is uniform for any value of the rates uab and uba- Interpreting A particles as 
particles and B particles as holes, we thus recover a known property of the ASEP [4, 5], 
namely that its stationary-state measure is uniform, irrespective of the bias. 

For three species {K = 3), we have /I3 = 5 and 5*3 = 3. There is indeed one single 
antisymmetric combination of pair numbers: 

Q = Nab- Nba = Nbc - Ncb = Nca - Nac- (4.6) 

There is accordingly a single condition on the six exchange rates for having a uniform 
stationary-state measure: 

uab + ubc + uca = uba + ucb + uac- (4.7) 

This condition is known in the context of the matrix-product formalism [16]. It can be 
checked that (4.7) is fulfilled by the rates of Table 8 in the absence of interactions 
{J a = Jb = Jc = 0). The only non-zero rates indeed read uab = wjj = 1/2, 
ucb = xij = (1 + q; + /3)/4, and uac = Vu = (1 — « — /5)/4, irrespective of /, J. 
In the CP-invariant case, one has uab = 1/2 and ucb = uac = 1/4- 
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Finally, for a large number of species {K ^1), asymmetric (driven, nonequilibrium) 
dynamics are far more constrained that symmetric (equilibrium) ones. Indeed the 
condition of having a uniform stationary-state measure roughly cuts off half the 
parameters, reducing their number from K{K—1) to Ax ~ -f^^/2, whereas for symmetric 
(equilibrium) dynamics the K{K — l)/2 rates are not constrained. The expression (4.5) 
shows that the difference A^ — Sk ^ K <^ Sk is relatively negligible for a large number 
of species. In other words, for the uniform stationary-state measure, the full space of 
nonequilibrium dynamics is hardly larger than the subspace of equilibrium dynamics. 

Example 2. 

Our second example still belongs to the realm of driven diffusive systems, albeit 
with multiple occupancies. The results below strengthen our conclusion and broaden 
its range of apphcability. 

Consider the class of dynamical urn models defined as follows. A^ particles are 
distributed among M sites around a ring, with multiple occupancies. Let A^^ be the 
number of particles at site m = 1, . . . , M. The system is subjected to the following 
stochastic dynamics. 

(i) a departure site d is chosen uniformly at random. 

(ii) a neighboring arrival site a is chosen as the right neighboring site {a = d + 1) with 
probability p, or the left neighboring site {a = d — 1) with probability q = 1 — p. 
(iii) a particle is transferred from site d to site a at a rate Wm which only depends on 
the occupancies k = N^ and I = Na oi the two sites. 

The relevant question in the present context is the following one. Under which 
conditions on the rates Wm is the stationary-state measure a product measure of the 
form 

P{C) = P{N,, ...,Nm) = ^^ PN, ■ ■ -PNm <5(iVi + ■ ■ ■ + Nm, N)1 (4.8) 

The answer to this question is known [17] (see [8] for a simple presentation). Consider 
first the case of an asymmetric dynamics {p ^ 1/2). The stationary-state measure is 
given by (4.8) if and only if the rates Wki obey the two conditions 

Pk+iPiWk+i^i = PkPi+iWi+i^k, (4.9) 

Wki - Wko = Wik - Wio. (4.10) 

The first condition (4.9) relates the rates Wki and the one-site factors pk of the 
stationary-state measure distribution. The meaning of this relation is clear: it just 
expresses P-related pairwise balance. The second condition (4.10), which does not 
involve the pk, is therefore more 'kinematic' than 'dynamical' in essence. 

The zero-range process (ZRP) corresponds to the particular case where the rates 
Wki = Uk only depend on the occupation of the departure site. The condition (4.10) 
is then automatically satisfied, whereas (4.9) yields the following relation between the 
rates Uk and the factors pk'- 

Uk = uj (4.11) 

Pk 
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for k > 1, where the constant uj fixes the time unit. 

The most general dynamical urn model with stationary-state product measure is 
hardly more general than the ZRP. Let us state the following result, skipping the proof. 
For a given product measure of the form (4.8), with prescribed factors pk, the general 
solution of (4.9) and (4.10) is entirely determined by the one-dimensional array of rates 
ttfe = Wko- One has indeed (with a^ = 0) 

1 ' 

Wkl = V Pk+mPl-m{<yk+m - O^l-m) ■ (4.12) 

PkPl m=0 

The rates ak are the rates at which an empty site (/ = 0) is refilled, by receiving one 
particle from a non-empty neighboring site containing k > 1 particles. 

In the case of a symmetric dynamics {p = 1/2), only the first condition (4.9) is 
requested [17, 8]. This relation expresses detailed balance. The resulting stationary 
state is therefore an equilibrium state. The condition (4.9) determines the rates W^i for 
k > I in terms of those for k < I. 

For a general dynamical urn model, the stationary product measure thus depends 
on the one-dimensional array of rates a^ in the asymmetric case, and on the two- 
dimensional array of rates Wki ioi 1 < k < I in the symmetric case. 

To sum up, the two above examples of driven diffusive systems corroborate the 
picture which emerges from the results of the present work. Asymmetric (driven, 
nonequilibrium) stochastic dynamics producing a given stationary-state measure are far 
more constrained (in terms of numbers of free parameters) than symmetric dynamics 
producing the same measure as an equilibrium measure. 
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